Linear Waves in Imperfect Relativistic MHD 

Fluid 



Sushil Kumar Singh and Subhash Kumar 
Department of Physics and Astrophysics, 
University of Delhi, Delhi-1 10007 (India) 



Abstract 

Linear relativistic anisotropic magnetohydro dynamic waves are investigated for an Imper- 
fect fluid. Consideration is carried out for the plasma with finite bulk viscosity. Conditions 
for growth and decay of these waves are obtained in special cases of isotropy and Chew- 
Goldberger-Low (CGL) state. 
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I. INTRODUCTION 



Relativistic plasma embedded in a strong magnetic field exhibits itself in a variety of interest- 
ing astrophysical events from pulsar winds to blackhole magnet ospheres QjJ El • Waves and 
instabilities for a single fluid with isotropic pressure have been investigated in the relativistic 
magnetohydrodynamic framework |UEj- The equation of state of plasma have also been subject 
of study when the strong magnetic field suppresses the equilibration of pressures parallel and 
perpendicular to itself (SJ. Plane magnetosonic waves have been studied for a single relativistic 
fluid with anisotropic pressures IB] . However, systems in which radiation can provide the 
mechanism for viscosity and heat conduction, the bulk shear can be effective. Fluid (imperfect 
fluids) systems with isotropic pressure and dissipative effects have been put to investigations 
by Wienberg 013]. 

The present paper intends to consider the effect of bulk viscosity on the plane wave propagation 
in relativistic magnetized plasma. The paper is organized as follows. In Sec. II we present the 
basic equations used in the present paper emphasizing the new energy-momentum tensor. In 
Sec. Ill the general dispersion relations are obtained and their derivation method is described 
briefly. In Sec. IV, to the first order of dissipation, the dispersion relation is discussed under 
various limits. 



We consider a plasma embedded in a magnetic field. Plasma fluid velocity and temperature 
both are allowed to be relativistic. Hereafter, we will use units where velocity of light c — 1. 
Let U a be the plasma 4-velocity and F a/3 be the electromagnetic field tensor. The plasma 
rest-frame electric and magnetic fields are defined [Zj as 



II. BASIC EQUATIONS 



E a =F af3 U, 



(1) 



B 



■ a 



(2) 
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respectively. The electromagnetic field tensor can then be written in the following form: 
where e al3lS is a completely antisymmetric tensor and e 0123 = 1. 

For a hydrodynamical system whose pressure, density and velocity vary appreciably over dis- 
tances of the order of mean free path, thermal equibrium is not strictly maintained and the 
fluid kinetic energy is dissipated as heat. The presence of weak space-time gradients in an 
imperfect fluid has the effect of modifying the energy- momentum tensor HJ • For a fluid system 
interacting with radiation and wherein bulk viscosity is the dominant dissipative component, 
the stress energy-momentum tensor is 

T ali = W x U a U p - W 2 r] aP - W z n a nP - (U^U 01 ^ - r/ Q/3 ), (4) 

where 

B 2 

W l = e+p ± + — , (5) 

W 2 = P± + ^, (6) 
B 2 

Wz=p±-p\\ + —. (7) 

Here B a = Bn a and the coefficient of bulk viscosity ( is positive definite. The plasma 4- velocity 
is normalized as U a U a = 1. The magnetic field unit vector n a obeys n a n a = — 1. The energy 
density e is in general a function of mass density p and the magnetic field strenght B, and can 
be put in the following functional form [Jj: 

e = pe(p,B) (8) 



The pressures are 



_ de de de . . 

p± = p ^ +pB dB = m- p -mB- e - (10) 

These equations are in their most general form and are representative of relativistic anisotropic 
state. 

Invoking the usual frozen-in condition E a = for the relativistic magnetohydrodynamic fluid, 
the electromagnetic dual tensor becomes 

G a ^ = B(n a U^ — n^U a ), (11) 

and the Maxwell equations yields 

G a l = 0. (12) 

The equations of motion of the fluid are contained in the particle and energy-momentum con- 
servation laws 0: 

J% = 0, (13) 
Tf = 0, (14) 

where J a (= pU a ) is the 4-current density. 

III. DISPERSION RELATION 



a 



The perturbed state is characterized by p — ► p+5p, B — > B+SB, U a — > U a +SU a , n a — > n a +<5n 
and plane wave solutions for the perturbations require that Sp, SB, 5U a , 5n a oc exp (ik a x a ). It 
should be noted that the condition U a U a = 1 and n a n a = — 1 yields the following constraints, 
U a SU a = n a Sn a = 0. 

It is assumed that the coefficient of bulk viscosity £ does not vary appreciably with variations 
in above system parameters (p, B, U a ). 

Substituting in Eqs.(12)-(14) d/dx a — > zA; Q , one obtains the following form: 
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kJJ a = 0, (15) 



kJG aP = 0, (16) 



k a 5T af3 = 0. (17) 

The frequency and parallel component of the wave vector of the plane waves are defined in the 
rest frame of the fluid as: 



u = k a U a , k\\ = -k a n a (18) 

respectively, and 

k a k a = s 2 = u 2 -k 2 = oj 2 -k\- k\. (19) 

The coefficients Wi in the expression for the energy-momentum tensor will be function of number 
density p and magnetic field strength B. Hence, 

SWi = —S P+ —SB. 

The four-vector defined as l a = e^^kpU^riv [7 j will be orthogonal to the independent vectors 
k a , U a and n a i.e., l a U a = l a n a = l a k a = 0. Multiplication of Eqs.(15)-(17) by l a , one obtains 



k\\ uj \f l a 5U a 
cuW 1 k\\W 3 J I l a 5n a 



= 0. (20) 
This yields the corresponding dispersion relation for the Alfven wave: 



oj = -r-ku, 



Wi 

(21) 



fP±-P\\ + €\ u 2 



Multiplying Eqs.(15)-(17) by the other three independent vectors k a , U a and n a gives the 
following matrix 
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where 



( u 














1 
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5lnp 







-k\\ 
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— UJ 












5lnB 







UJ 








—k\\ 


1 






k a Sn a 




At 


B x 





k\\W 3 





w x 






UJn a 




A 2 


B 2 


2k\\W 3 








^iCujWik 2 






n a 5U a 




\ A 3 


B 3 


w 3 





ujW x 


-iCk\\ 


) 


\ 


k a 5U a 


) 



A, 



d 



dlnp 



2), 



(22) 



a 3 = h-^(w 2 -W 3 ) 



d 

A 2 = ujA x + k\\A 3 + k\-—W 2 

dlnp 



B^uj^( Wi -W 2 ) 



B 3 = ki 



d 



dlnB 



(W 2 -W 3 ), 



B 2 = ujB 1 + k\\B 3 + k 







dlnB 



-Wo. 



The corresponding dispersion relation for the slow and fast magnetosonic modes are given by 



■ 2 d(W 1 -W 2 ) u2 d(W 2 -W 3 ) 



dlnp 



+k- 



dW 2 
dlnp 



k?AW 3 + 



dlnp 
d(W 2 - W 3 ) 



k 2 W 3 - uj 2 W x + k 



dW 2 



dlnB 



+ uj 2 [W 1 



dlnB 
d{W 1 - W 2 
dlnB 



-i(uj 



ki kiW 3 - uj'W x )+kUuj 



,d{W 1 -W 2 ) 
dlnB 



11 dlnB ) 



0. 



(23) 



In the limit of bulk viscosity £ — > (perfect RAM fluid), this equation reduces to the dispersion 
relation obtained by Gedalin [7j. 
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IV. DISCUSSION 



Unmagnetized Plasma : 

In the absence of magnetic field, the energy density e = e(p) and isotropy prevails i.e., W 3 = 0; 
and to the first order of dissipative term , we obtain 

uo = uj + iLk 2 , (24) 

where the decay characteristic length 

c 



2(e + pY 

and the frequency of the sound waves with the bulk shear being neglected 



(25) 



uj = k 



dW 2 
dlnp 



2 



The amplitude of the sonic waves decays at the rate T = Lk 2 . 
Magnetized Plasma : 

(i) In the special case of propagation of magnetosonic modes in the direction of the magnetic 
field, the frequency is found to be 

u = uj + iL\\k 2 , (27) 



with the decay characteristic length 



L \\ = -T- 1 7, (28) 



and the viscosity neglected frequency 

Uq = k 



d(W 2 -W 3 ) i 



dlnp 



dlnp 



2 



(29) 
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It is to be noted that L\\ is always positive. This implies that the wave decays for all possible 
values of the system parameters. 

(ii) For the magnetosonic waves propagating perpendicular to the ambient magnetic field, the 
frequency is 

lj = ujq — iL±k 2 , (30) 

where the spatial scale of decay 



2 \ dlnB Z 4tt > 



T _ S V dlnB ~ 4tt / /oi\ 
-L ~ o/_ , , B*\ 9e > 



and the frequency in the absence of dissipation 



a(H/i-VF 2 ) aw 2 I dW2lyy _ d(W 1 -W 2 ) -\ N i 

dlnp dlnB dlnp L ^ dlnB 



d(W 1 -W 2 ) 

8p± 



2 



cfeB ^ 4tt ^ 3p- ^ dlnB > \ 2 

olnp 1 



= * f^— gs • (32) 

In this case we observe that the magnetosonic modes can grow or decay accordingly as Lj_ > 
or L± < 0. In the latter case the system tends to return to its equilibrium state while in the 
former instability grows. A state characterized by L± = represents propagation of harmonics 
in the direction transverse to the magnetic field unlike the unmagnetized-isotropic state wherein 
the harmonics can propagate only upto its few characteristic lengths. 

Shrauner JU] proposed a generalized polytrope model for the MHD set of equations for con- 
sidering the collisional, collisionless and the transitional regimes of plasma state. According to 
Shrauner, the anisotropic pressures can be written as 

P]l =C llP »B- a , P± = C ±P 5 B^ (33) 

where C\\ and C± are positive constants. The polytropic indices u, a, 5 and (3 are positive 
constants and these polytropic relations can be reduced in certain special cases to well known 
equations of state for the pressure. From equations (8), (9), (10) and (33) we get the generalized 
energy density 

e = aip + a 2 p\\ + a 3 p_L, (34) 
8 



where ai, a 2 and 03 are positive constants. The condition for decay of the transverse modes i.e., 
L±_ < translate into 

B 2 

a 2 ap\\ > a 3 (3p± + — . (35) 



For ai = 1, a 2 = 1/2 and 03 = 1 we recover the energy density for a plasma state characterized 
by nonrelativistic temperatures; e = p + \p\\ + p±. Further, if v — 3, a — 2, S — 1 and (3 = 1, 
then the system represents the well known CGL ^T] state with py = C\\p 3 /B 2 and pj_ = C±pB. 
It is evident that the condition for decay transverse modes in the CGL plasma reads as 

B 2 

P\\>P± + —- (36) 
This is the usual criterion for firehose instability in MHD. 

We would like to point out that a simple gas of structureless point particles will have negligible 
bulk viscosity in the extreme-relativistic or nonrelativistic limits. However, this need not be 
the case for a fluid composed of a mixture of highly relativistic and nonrelativistic particles. 
It is well known that the exchange of energy between translational and rotational degrees of 
freedom gives ordinary diatomic gases an appreciable bulk viscosity. Although the present 
analysis assumes that coefficient of bulk viscosity ( does not vary appreciably, it does depend 
on the dynamical variables of the system in the general form ( = 4&T 4 t[| — (^) p ] 2 . Here T is 
plasma temperature, b is of the order of Stefan-Boltzmann constant and r is the free mean time 
of radiation quanta [9 . In view of this, the present work is an indicative analysis of plasma as 
an imperfect fluid and can be explored further. 
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